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Abstract. We present magnetic properties of the three-band Hubbard model in the para- and antiferro- 
magnetic phase on a hypercubic lattice calculated with the Dynamical Mean-Field Theory (DMFT). To 
allow for solutions with broken spin-symmetry we extended the approach to lattices with AB-like structure. 
Above a critical sublattice magnetization md ~ 0.5 one can observe rich structures in the spectral-functions 
similar to the t-J model which can be related to the well known bound states for one hole in the Neel- 
background. In addition to the one-particle properties we discuss the static spin-susceptiblity in the para- 
magnetic state at the points q = and q = (ji,-n,Ti, ■ ■ •) for different dopings 8. The 5-T-phase-diagram 
exhibits an enhanced stability of the antiferromagnetic state for electron-doped systems in comparison to 
hole-doped. This asymmetry in the phase diagram is in qualitative agreement with experiments for high-T c 
materials. 

PACS. 71.27.+a Strongly correlated electron systems - 71.30.+h Metal-insulator transitions and other 
electronic transitions and 75.10.-bGeneral theory and models of magnetic ordering 
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1 Introduction and Model 

One of the few undisputed facts about high-T c materi- 
als is that all undoped high-T c compounds are insulators 
with antiferromagnetic ordering in the Cu02-planes at low 

. enough temperatures jjj . Doping of the systems leads to 

' a strong suppression of the antiferromagnetic order and 
eventually superconductivity sets in. While the explana- 
tion of this transition and the proximity of antiferromag- 
netism and superconductivity surely is the most fascinat- 
ing aspect in the high-T c compounds, there are a variety 

' of other peculiarities that call for an explanation. One of 
these side aspects is the observation that, although the 
basic scenario is the same, hole and electron doped mate- 
rials show an apparent qualitative difference in that the 
magnetic phase in the latter appears to be much more 
stable for the latter (cf. Fig. |l|). 

The physics in the insulating phase is well described 

1 by a Heisenberg model since charge fluctuations are sup- 
pressed by strong local correlations and, to a first ap- 
proximation, one is left with the Cu-spin degrees of free- 
dom only. The exchange parameter of the resulting effec- 
tive Heisenberg model is obtained from standard superex- 
change processes ||||, induced by virtual hopping of a 
hole from one Cu-ion to the neighbouring one over a non- 
magnetic O-ion. 

At any finite doping one has to consider at least also 
the charge degrees of freedom on the copper sites and 
would then be left with the usual one-band Hubbard or 
t-J model to describe the interplay between magnetic ex- 
change and itinerancy. However, this scenario completely 
neglects the existence of the oxygen sites. That they are 
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Fig. 1. Experimental phase diagram for high-T c compounds 
(taken from 111]). 



indeed important, at least for the magnetic properties, 
can be seen from the following qualitative argument: Due 
to the strong local correlations at the Cu-sites, additional 
doped holes mainly occupy O-sites. The spin of the hole 
at the O-site induces an effective ferromagnetic interac- 
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tion j4j between the neighbouring Cu-spins, so that the 
antifcrromagnetic ordering is strongly suppressed with in- 
creasing hole-doping. In the case of electron doping, on the 
other hand, the additional particle has to go to the copper 
sites due to Pauli's principle, which effectively means that 
free spins are removed from the system. This obviously 
also leads to a suppression of magnetic order, but in a 
weaker fashion. Thus, in order to obtain a more realistic 
description of the physics of high-Tp compounds one has 
to take into account also the oxygen degrees of freedom 
and therefore one-band models like the standard Hubbard 
or t-J model become inadequate. 

The simpliest model which includes both the effect of 
strong correlations and the influence of the oxygen sites, is 
the three-band Hubbard model or Emery model ||. The 
three-band Hubbard Hamiltonian reads 

H = £ d d L d ™ + £ pP]aPj° 

+ Y *« ( d \*pj° + h - c ) +J2 Ud AA > (!) 

where the vacuum is denned as all orbitals in ([[]) filled with 



electrons 
a Cu3d„2 



With this convention d' i(T (plj a .) creates a hole in 
_ y 2-(0 2p x / y )-orbital at site with spin a and 
Ed{£ P ) are the corresponding on-site energies. Uj = ±f de- 
notes the nearest-neighbour hopping matrix- element be- 
tween Cu- and O-sites and Ud stands for the Coulomb- 
interaction of two holes, residing at the same Cu-site i 
with number operators nf a . 

In this paper we want to study magnetic properties of 
the Hamiltonian ([!]) in the framework of the DMFT. In 
this theory the dynamical renormalizations of the one par- 
ticle properties become purely local Hff) , so that they can 
be obtained from an effective impurity problem coupled to 
a self-consistent medium. Due to the additional orbital de- 
greee of freedom in (Q) the mapping on the corresponding 
effective impurity model is not unique. In order to treat 
local spin and charge fluctuations between the Cu-site and 
the surrounding O-sites better than on a mean-field level 
we use the approach, developed in ref. ||, where a cluster 
of one Cu d-orbital and a normalized bonding combination 
of the four surrounding O-orbitals is coupled to the effec- 
tive medium. This method indeed leads to the anticipated 
physics on the one-particle level, namely the formation 
of a low-lying singlet state - the Zhang-Rice single t O ] 
and at half filling to a charge-transfer insulator |p|,p|JlC[f in 
contrast to the Mott-Hubbard scenario for the one-band 
model. So far, however, only the paramagnetic state has 
been studied in ref. || . In ref. the metal- insulator (MI) 
transiton was studied with Quantum Monte Carlo in the 
context of the DMFT. There the model was also embed- 
ded on a bipartite lattice in order to take into account the 
antiferromagnetic symmetry breaking, but in this article 
the attention was called to the Mi-transition. 

In order to obtain the phase diagram or look at the 
behaviour in the antiferromagnctically ordered state the 
method has to be extended to allow for the calculation of 
susceptibilities or solve the DMFT equations for lattices 



with AB-like structure, respectively. A short review of this 
generalization together with the technique of calculating 
the magnetic susceptibility will be given in the next sec- 
tion, followed by the discussion of our results in section 3. 
The paper will conclude with a summary and outlook in 
section 4. 



2 Method 

2.1 The DMFT for the three-band model 

Let us begin by summarizing the basic concepts intro- 
duced in H for the three-band Hubbard model. In order 
to construct the DMFT for a Cu-0 plaquette, it is con- 
venient to introduce the Fourier transform of the kinetic 
part of the Hamiltonian (l]) after generalization to d di- 
mensions, which then reads Iql 



H = h a (k) + u d « + H n 



(2) 



with 



Mk) 



Sdd^dk, 



2dtj k (d^p^ + h.c 



(3) 



Here, c?ko- is the Fourier transform of di a and pkc is the or- 
thonormalized Fourier transform of the hybridizing com- 
bination of the oxygen orbitals surrounding a given cop- 
per site |U|]- The d — 1 linear combinations, which are 
orthogonal to pko- were collected into -ff n on-bond and will 
be dropped in the following, because they are decoupled 

from the remainder of the system. Finally, 7k is given by 

d 

% ' ' 
rescaling 



1 — 4 2 cosA;„. In ref- S it was shown, that the 



with 7 k = w 1 



leads to a nontrivial limit for d 
tion in the DMFT now takes the form 



2^ 7k -> 2t* 7 £ (4) 

const. = 1 andek = Y^t=i cos kv 
oo. The d-Green's func- 



Gt(z) = 



z-e d ~ E*{z) ~ 



At* 



V2d K 



(5) 



The new ansatz by Schmalian et al. was to write the local 
d-Green's function to be of the form H 



k 

~z-e d - S d a {z) ~ 



(6) 



At* 



In the DMFT the effective Cu-0 cluster lives in a so-called 
effective medium, defined via 

At* 2 



g a {z)- 1 =G d a {z)- l +£ d a {z) = z- 



z — £ p — A a (z) 



(7) 
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Note that in this form the coupling to the rest of the sys- 
tem, which is described by A a (z) happens through the p- 
states only This representation of the local Green's func- 
tion is obviously not unique. One could also choose a rep- 
resentation for the local e?-Green's function of the form 

r i ~ 1 

G^(z) = z — £d — S^(z) — A a (z) where the resonance 

at z = e p is included in A G (z). But from a numerical 
point of view the form (^|) is more convenient because the 
singularity at z = e p is not included in the hybridization 
function which therefore becomes smooth as a function of 
frequency. The form (^J) of the local Green's function is 
just the Dyson equation of an effective impurity problem 
consisting of one d- and one p-orbital, where only the p- 
orbital hybridizes with the conduction electrons (see eq. 
13 in ref. §). 



which are different in the antiferromagnetic state. Fur- 



thermore 77^ = ^ ^7*. ± 7k+cjJ ■ For the <i-components 
of the Green's function matrix we finally obtain 

/ 

g°(k,z) = C ka 



(12) 



with Cka 



/ 2 
_ ( 4t* £ k 

AK '- r \V2dH P 
it* 



£p = Z + fj, 



e p and 



= £ d — 2Ja/b~ "7 — ■ The l° ca l Green's function 



J A/B 



is obtained by taking one of the diagonal elements and 
summing over k with the result 



2.2 DM FT for the Neel state 

In the antiferromagnetic phase the period of the unit cell 
of the lattice is doubled due to the reduced translational 
symmetry. Consequently, the volume of the magnetic Bril- 
louin zone (MBZ) is reduced to one-half of the volume in 
the paramagnetic state and the vector Q = it, 7r, • • •) 
becomes a reciprocal lattice vector. These changes in the 
symmetries of the system can be simply taken into account 
by introduction of an AB-sublattice structure and re- 
formulating the theory on an enlarged unit cell containing 
exactly one A- and one B-site. Since this procedure does 
not affect the local two-particle interaction in the Hamil- 
tonian (^) we will concentrate on the kinetic part for the 
derivation of the resulting Hamilton matrix. We first split 
the kinetic part of (||) in the following way: 



H= (Mk) + Mk + Q)} 



(8) 



keMBZ.o- 



Note that the k-sum runs over k-points in the reduced 
Brillouin zone only! Rewriting (||) in terms of the linear 
combinations 



(9) 



PA/Bka = ^75 (Pkcr ±Pk+Q<j) , 

acting on the A- or B-sublattice, respectively, one obtains 
H= d a (k)^ka • (10) 

keMBZ.o- 

For simplicity we introduced a spinor notation for the op- 
erators = (<4 kff pi kCT 4jkaPBk CT ) and thc Hamilton 
matrix on the sublattices 



yo tj+ 



IL(k) 



V 















77 k £ d 




J7 <T 77+ 
B u k 



(ii) 



'k u xx k £ P I 

The quantities denote the local self-energies due 

to the two particle term in (||) on A/B-sublattice sites, 



Gd,A/B{ z ) 



ds Po(e)- 



B/A 



(13) 



where the density of states p (e) corresponding to the 
dispersion £ k was introduced. In the paramagnetic state, 
where £1 = one immediately recovers the result of 



section 2.1 



In the antiferromagnetic state it is sufficient to per- 
form the calculations for the A-sublattice only due to the 
additional symmetry jl2| G% A — G^ B and use the spin- 
index for book-keeping. The actual calculation is now a 
straightforward extension of the method used in ref. @] 
for the paramagnetic phase. The local nature of the self- 
energies allows the mapping of the lattice problem on an 
effective impurity-problem, consisting of a <i-orbital and 
the orthonormalized hybridizing combination of the four 
surrounding p-orbitals, coupled to the effective medium, 
which is described by the propagator (Q) and has to be de- 
termined selfconsistently. Again, the coupling to the sur- 
rounding clusters is assumed to happen through the p- 
states only. 

The remaining local problem is solved with the resol- 
vent method p3UlJ] and an extended version of the so 
called Non Crossing Approximation (NC A) [Q , where the 
16 local eigenstates of the impurity are coupled through 
the hybridization- function A cr (z) J|. 



2.3 On the calculation of the magnetic susceptibility 

On the one-particle level one can obtain magnetic proper- 
ties by applying a staggered magnetic field and calculat- 
ing the sublattice magnetization. This technique is very 
tedious so that we used another method for calculating 
the magnetic phase diagram. 

In addition to the one-particle properties the DMFT 
also allows to calculate two-particle correlation functions, 
e.g. the magnetic susceptibility consistently. In analogy to 
the one-particle case the two-particle self energy becomes 
purely local in the limit d — ► 00 JlJ,[l5| ■ This enables us to 
extract the two-particle self-energy from the effective local 
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problem p5|jlq , |17| and use it to determine the two-particle 
correlation function for the lattice. Since the local two- 
particle propagator is a function of three frequencies in the 
most general case the algorithm works best for Matsubara 
frequencies, because all quantities can be represented as 
matrices in this case. For details of the method see e.g. 
ref. @. 

The choice of the cluster as effective impurity and the 
finite value of Ud results in 16 local eigenstates. This leads 
to a huge number of diagrams for the local two-particle 
propagator, which have to be calculated as functions of 
three frequencies and summed up numerically. Although 
the problem of generating the correct diagrams for the 
local two-particle propagator can be automated and han- 
dled by the computer, the remaining numerical task is still 
formidable and restrict our calculations to the evaluation 
of the static susceptibility for the time being. Nevertheless, 
the study of the dynamical susceptibility is in principle 
also possible |Q and will be the subject of a forthcoming 
publication. 



3 Results 

3.1 Susceptibility and phase diagram 

Let us start the discussion of our results with the magnetic 
susceptibility. We calculated the static magnetic suscep- 
tibility of the Cu-spins in the paramagnetic phase at the 
points q = and q = Q, which give the homogeneous 
and staggered susceptibilities, respectively. For the param- 
eters of the three-band Hubbard model we have chosen 
Ud = 2 A = 7t* , where the charge transfer gap A is de- 
fined by A — e p ~Ed- Fig. ||a shows typical results for these 
two susceptibilities as a function of temperature T for a 
hole-doping 5 = rid + n p — 1 = 0.1. For the above choice 
of Ud and A we observe a finite and slowly varying ferro- 
magnetic susceptibility Xf(T, S), which does not show any 
tendency towards an instability in the calculated region of 
temperatures and dopings. The antiferromagnetic suscep- 
tibility Xaf(T, S), on the other hand, varies strongly as a 
function of temperature and diverges at a finite tempera- 
ture T — Tm- Fig. [2)d shows the inverse staggered suscep- 
tibility for the same parameters. As expected we find the 
linear variation of Xaf' which is typical for a mean-field 
theory. By calculating the inverse susceptibility for differ- 
ent dopings 5 we obtain the (5-T-phase-diagram, shown in 
Fig.pi. Note, that half-filling (<5 = 0) does not coincide with 
the Mi-transition. This is so because for the used param- 
eter values of Ud and A the metal-insulator-transition is 
shifted towards larger hole-filling values n > 1. Obviously 
the highest value for the Neel temperature is achieved at 
the mctal-insulator-transition and not at half-filling. As 
already observed for the one band model [ pj|Jl^1 the anti- 
ferromagnetic phase is strongly suppressed upon doping. 
However, in contrast to the former case one recognizes 
a pronounced asymmetry in the Neel temperature with 
respect to hole- and electron-doping. This stronger sensis- 
tivity of the antiferromagnetic ordered state in the case of 
hole doping compared to electron doping is qualitatively 
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in good agreement with experiments (see Fig. [j]) and is a 
direct consequence of the oxygen degrees of freedom. The 
present large D treatment however neglects short range 
order phenomena so that the argumentation of ref. [|4| 
concerning the frustration of the AF-order by hole doping 
does not hold on this level. But nevertheless the oxygen 
degrees of freedom cause a particle hole asymmetry in the 
half filled three band model which yields the calculated 
asymmetry in the phase diagram in this large D approach. 

In our calculations the ordered phase is stable up to 
S 0.18 for low enough temperatures (cf. Fig. ||). Ex- 
periments and other theoretical calculations show much 
stronger suppression of the antiferromagnetic ordering with 
doping Q. The tendency to overestimate the magnetic 
phase boundary is typical for mean field theories, since 
they completely neglect fluctuations, which strongly renor- 
malize transition temperatures. 



3.2 Spectral functions in the ordered phase 

With the generalized equation (|l^) it is also possible to 
perform calculations in the antiferromagnetic phase. To 
allow for solutions with finite sublattice magnetization 



A, (co) 



rrid 



A/B A/B 



'<i\ 



we apply a small symmetry-brea- 
king staggered magnetic field hlji) = he lCi ' ri in z-direction 
at the beginning of our iteration procedure, which is turned 
off after a few iterations. 

In the following we concentrate on the <i-part of the 
spectrum A%(lj) = — — ImG^ A( u>Jr ^) on the A-sublattice, 
since the p-part shows exactly the same features as the d- 
part, only with different spectral weights for the various 
bands. Fig. f|a,b,c shows the typical behaviour as the tem- 
perature T is lowered for fixed values of the parameters 
U d = 2Z\ = It* at finite doping 6 = 0.015. 

For j3 = 18/t* the system is in the paramagnetic state 
(cf. Fig. [Ja], where the general equation ( |l3|) reduces to 
the form (||). We find the same result for the <i-part of 
the spectrum as in ref. [^J, where a detailed discussion 
of the various bands concerning their doping dependence, 
transfer of spectral weight and the evolution of coherent 
quasiparticles near the Fermi-energy can be found. Let 
us just briefly mention the important low energy parts, 
namely the so-called lower Hubbard band at to < in the 
spectrum of Fig. [|, which has mainly d-character, and the 
Zhang-Rice band right above the gap, which is generated 
by the singlet combination of the p-and d-states on one 
plaquette B. Note that S > although the chemical po- 
tential is located in the lower Hubbard band. Therefore 
the Mi-transition oc curs at larger filling values as already 
mentioned in section I 



3.1 



With decreasing temperature the 
system enters the antiferromagnetic phase and the spec- 
tral functions of up and down spin become inequivalent 
(cf. Fig. ^b), yielding a finite sublattice magnetization nid- 
Note that the major effect is a transfer of spectral weight 
from the minority spin to the majority spin. In addition 
the peaks in the spectra are slightly shifted in energy with 
respect to each other. This effect can be ascribed to an in- 
ternal molecular field, generated by the finite sublattice 




10.0 



Fig. 4. d-spectral function for the majority (full curves) and 
minority spin (dashed curves) for Ud = 2/1 = It*, fixed dop- 
ing 5 — 0.015 (below the Mi-transition (cf. Fig|3|)) for various 
temperatures (3 = 18/t* (a), /3 = 20/t* (b) and f3 = 22/i* (c) 
and for 5 = 0.063 (above the Mi-transition) and /3 = 22/t* (d) 



magnetization. Therefore, measuring the energy shift one 
can calculate the internal molecular field and from this 
the exchange parameter J of a corresponding tJ-model. 
For still lower temperatures the sublattice magnetization 
increases and above a value of md ~ 0.5 a pronounced 
multipeak structure is evolving (see Fig. ^c). Fig. ||d shows 
the spectral function for the same system parameters and 
sublattice magnetization as in Fig. |jc but with the chem- 
ical potential right above the gap. Also in this regime the 
same multipeak structure occurs and the spectral func- 
tion shows little difference to Fig. |]c. Only the peaks next 
to the chemical potential have more spectral weight com- 
pared with Fig. I]c. 

These regular resonances were previously found in DMFT 
calculations of the t J-model |n| . There the multiple peaks 
could be related to bound states of one single hole in the 
Neel background. In ref. |^0| it was shown, that this spe- 
cial problem can be solved exactly for d = oo and T = 
within the tJ-model. The most important physical aspect 
is, that the moving hole feels a binding potential propor- 
tional to J, growing linearly with the distance from its 
starting point due to the breaking of antiferromagnetic 
bonds during its motion J|(]Jl9| . This linear potential leads 
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Fig. 5. Comparison of the low energy part of the d-spectrum 
with the exact results (vertical lines), given by equation 



to a sequence of discrete poles at frequencies 

a n t 



-It 



J* 

Y 



(14) 



as spectrum for the one particle excitations. Here, the a n 
denote the zeros of the Airy function Ai [At/ J*), and the 

renormalized parameters i and J* are given by i = t\/2d, 
J* = J2d. These exact results can be compared directly 
to the resonances, found in the DMFT calculations for the 
tJ-model Since the model used in our calculations is 
fundamentally different from the t J-model and also from 
the one-band Hubbard model, the relevance of this physi- 
cally intuitive picture to it and especially the proper choice 
of the parameters for an effective t J-model to describe the 
low energy properties is not clear a priori. The approach 
chosen here is to fix the hopping to i 2 /A, which reproduces 
the free bandwidth. In our case we determine an effective 
exchange interaction J* from the energy shift AE — J*m ( i 
of the bands. Note that this is just the energy shift of the 
spin-up and spin-down bands of a corresponding t J-model, 
treated on a mean-field level j?],[l9|j2l[] . Another possibilty 
to obtain the exchange integral J is to use the result of 
a Schrieffer-Wolf transformation of the 3-band Hubbard 
model, see e.g. However this transformation holds only 
for large values of Ud and A, so that we do not expect this 
procedure to give a meaningful result for our parameter 
values. 

Fixing the paramters in eq. ( |l4| ) as discussed above, 
we can indeed directly compare our results with the dis- 
crete spectrum ([jj]) . Fig. || shows some examples for the 
fit of the low energy part of the <i-spectrum A^(uj) by 
the discrete spectrum ( |l4|) at fixed doping 5 — 0.015 and 
sublattice magnetization rrid = 0.60 for various parame- 
ters Ud — 2A. Note, that the energy scales of up- and 
down-spin in Fig. || are already shifted by ±AE/2 respec- 
tively, so that the resonances of majority- and minority- 
spin bands coincide. We find quite good agreement with 
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minority spin (dashed curves) for Ud = 2A — It* and (3 = 
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the distance of the peak positions. The broadening is ex- 
pected to result from finite temperature, sublattice mag- 
netization and doping effects |L9|. This means that the 
tJ-model with proper choice of the parameters t and J 
seems to reproduce the low energy one-particle dynam- 
ics of the three-band Hubbard model in d = oo correctly, 
even in the antiferromagnetic state. In addition, the basic 
physical picture for the multipeak structures observed for 
low temperatures appears to be the same as in the simple 
one-band models. 

In order to gain more insight in the effect of doping 
on these multipeak structure we investigated the spec- 
trum at larger doping far away from the Mi-transition. 
Fig. | shows the results for the d-part of the spectrum 
for the same system parameters and sublattice magneti- 
zation as in Fig. |Jc,d at (3 = 50/t* but at larger doping 
5 = -0.08 (a) and 5 = 0.13 (b). In the electron (Fig. 
^a) as well as in the hole doped regime (Fig ||b) only the 
resonances next to the chemical potential survive. Due to 
the larger doping there are more electrons/holes in the 
system whose paths can intersect and restore the anti- 
ferromagnetic background. Therefore the electrons/holes 
become more mobile and the resonances at higher energies 
are washed out. 

In finite dimensions the string picture for one hole in 
the antiferromagnetic background no longer holds and is 
correct only up to order jp Q due to the possibility of 
paths which intersect and touch themselves ^2|. Second, 
fluctuations become more important which can restore the 
antiferromagnetic background. Thus in low dimensions we 
expect that the multipeak structure at finite doping will 
dissappear. 



4 Summary 

In this paper we presented results for the magnetic proper- 
ties of the three-band Hubbard model in the limit of high 
spacial dimensions. These were obtained in the framework 
of the Dynamical Mean Field Theory, which enabled us 
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to calculate the one particle spectrum as well as two par- 
ticle correlation functions, namely the magnetic suscep- 
tibility. From this we evaluated the 5-T-phase diagram, 
which shows strong suppression of the antiferromagnetic 
state upon doping. In contrast to one-band models the 
ordered state is found to be more sensitive upon doping 
in the case of hole doping in comparison to electron dop- 
ing. This asymmetric behaviour is qualitatively in good 
agreement with experiments. The spectral function for 
single particle excitations in the antiferromagnetic phase 
shows pronounced features above a sublattice magnetiza- 
tion ma ~ 0.5. These structures are similar to those found 
in the t J-model for the special case of one single hole, mov- 
ing in the Neel background and can be understood by the 
binding of one hole in a string potential. A quantitive fit 
of the spectral functions by the exact results for the spe- 
cial case for the t J- model shows quite good agreement, so 
that the tJ-modcl seems to reproduce the correct physics 
of the three-band Hubbard model as long as one is only 
interrested in the low energy one-particle physics. 

In low dimensions fluctuations become more important 
which will destroy the multiple peaks found in the spectral 
function for d = oo. Thus these peaks have not yet been 
observed in experiments. 
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